Kondo physics in transport through a quantum dot with Luttinger hquid leads 
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We study the gate voltage dependence of the linear conductance through a quantum dot coupled 
to one- dimensional leads. For interacting dot electrons but noninteracting leads Kondo physics 
implies broad plateau-like resonances. In the opposite case Luttinger liquid behavior leads to sharp 
resonances. In the presence of Kondo as well as Luttinger liquid physics and for experimentally 
relevant parameters, we find a line shape that resembles the one of the Kondo case. 



Electron correlations can strongly alter the low-energy 
physics of many-electron systems. The Kondo effect and 
Luttinger-liquid (LL) behavior are two of the most promi- 
nent examples, both affecting electron transport through 
a quantum dot coupled to one-dimensional (Id) leads. 

For noninteracting leads the appearance of Kondo 
physics was investigated theoretically for the (two-lead) 
single impurity Anderson model (SIAM).P] At small 
temperatures T and for sufficiently large U/ A the Kondo 
effect leads to a resonance in the linear conductance 
G(Vg) with an unusual broad plateau-like line shape re- 
placing the Lorentzian resonance (of width 2A) known 
from tunneling at J7 = 0. Here U denotes the local in- 
teraction on the dot, A = -I- A/j measures the hy- 
bridization of the dot and the (left and right) lead states, 
and Vg is the gate voltage applied to the dot region. On 
resonance the number of dot electrons is odd implying 
a local spin-1/2 degree of freedom on the dot that is 
responsible for the Kondo effect. For T — > the res- 
onance height approaches 2(e^//i) 4AiAfl;/(Ai -I- A^)^, 
i.e. the unitary limit for symmetric dot-lead couplings, 
and its width is of order u'^,^ For the SIAM at T = 0, 
G is proportional to the one-particle spectral weight of 
the dot at the chemical potential ii. 0] Varying Vg within 
an energy range of order U the Kondo resonance of the 
spectral function is pinned at (close to) fi and has a fixed 
height which explains the broad plateau- like resonance 
in G(Vg). A series of transport experiments on quantum 
dots was interpreted in the light of these results. 5] 

The line shape of G{Vg) is equally strongly affected 
by the correlations in the Id leads if the Kondo effect is 
suppressed. This can be achieved considering one of the 
three cases: no spin degeneracy on the dot; no interaction 
on the dot; spinless fermions. The low-energy physics 
of Id wires of interacting electrons is characterized by 
a vanishing quasi-particle weight and power-law scaling 
of correlation functions known as LL behavior. :6:] In the 
case of spin-rotation invariant interactions (and spinless 
fermions) all exponents can be expressed in terms of a 
single LL parameter Kp that depends on the interaction, 
the filling factor n, and other details of the model con- 
sidered. For repulsive interactions Kp < 1. At T = 
and for finite LL leads the resonances at Vg have ap- 
proximately Lorentzian shape. For a spinful model and 



symmetric dot-lead couplings G(VJ) = 2e^ /h. At asymp- 
totically large N, where N is the length of the LL leads, 
the width w vanishes as a power law, in strong contrast 
to the broad resonances induced by the Kondo effect. 
At asymmetric couplings there are still resonances but 
limjv^oo G = also at VJ. 7] 

The problem of a single spin-1/2 coupled to a LL was 
investigated generalizing the Kondo model. However, 
the very interesting question of the resonance line shape 
resulting from the competition between the two correla- 
tion effects was not addressed so far. Here we will inves- 
tigate this fundamental issue. We study three cases using 
an approximate method that is based on the functional 
renormalization group (fRG): 9] (a) noninteracting leads, 
interacting dot; (h) LL leads, noninteracting dot; (c) LL 
leads, interacting dot. We focus on T = 0. Unless oth- 
erwise stated we consider symmetric dot-lead couplings. 
For case (a) we reproduce the pinning of spectral weight 
at fi and thus the plateau- like resonance. For case (b) 
we confirm the expected LL line shape of the resonances. 
Some emphasis is put on the two-particle backscatter- 
ing that for spinful particles plays an important role on 
intermediate length scales and was not investigated so 
far. The results show that the aspects of Kondo and LL 
physics essential to answer the above question are cap- 
tured by our method. For case (c) we show that for LL 
leads of experimentally accessible length the line shape 
resembles the one of case (a). Our results indicate that 
the plateaus vanish for — *■ oo. 

The model we investigate is given by the Hamiltonian 

H = i7kin + ^^int + ^^bar + -ffgatc , (1) 

with the kinetic energy 

oo 

o-=T,i j = -oo 

the interaction 

N N-1 
i?int = E^j"j.T%i + E E ^i"j.'^"j + b<T ' 

where fij^^ = Uj^^ — v (see below), the tunnel barriers 

i?bar = (t-0 E (d-l.^S,.+d.-S. + l,.+H-C.) , 



2 



and the part containing the gate voltage 
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Standard second quantized notation is used. The inter- 
action is restricted to a finite part of the wire {j G [1, A'']) 
corresponding to an experimental setup where the LL 
wires are connected to (higher dimensional) Fermi liquid 
(FL) leads. The model with 7?kin + ^^int and interaction 
on all sites is known as the extended Hubbard model. 
Away from half- filling, that is for n ^ 1, it is a LL at least 
for sufficiently weak repulsive interactions. Very accurate 
results for Kp were recently obtained numerically. 0| We 
allow the onsitc Uj > and nearest-neighbor Uj > 
interactions to depend on position. By turning on the 
interaction adiabatically over a few 10 lattice sites we 
can suppress any electron backscattering at the FL-LL 
contacts, which we are here not interested in.^^ The 
constant bulk values of Uj and Uj are denoted by U 
and U' , respectively. We can choose different interac- 
tions in the leads and on the dot sites (j S [jhjr] with 
Nd — jr — ji + 1) and thus model the cases (a) to (c). 
As long as ji and jr are sufficiently far away from the 
contacts at sites 1 and N the position of the dot does 
not play a role. In TJjnt we shifted rij.o- by the parameter 
h' = v{n, U, U') which is chosen such that on sites 1 to 
iV, but excluding the dot region, the average density n 
acquires the desired value, (llj This is important as Kp 
depends on n. The tunnel barriers are modeled by re- 
duced hoppings t' < t across the bonds linking the sites 
ji - 1, ji and jr,jr + 1- 

At temperature T = the linear conductance of the 
system described by Eq. |^ can be written asjl2j 
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with the effective transmission \t{e,Vg, N)\'^ ~ (4t^ — 
[e + fi]'^)\Gi,N{£,Vg, N)\'^ . The (spin independent) one- 
particle Green function G has to be computed in the 
presence of interaction and in contrast to the noninter- 
acting case acquires an N dependence. In the notation 
used in the following we suppress the argument N in G. 
To determine Q we use a recently developed fRG scheme. 
Staring point is an exact hierarchy of differential flow 
equations for the self-energy matrix and higher or- 
der vertex functions, where A G (oo,0] denotes an in- 
frared energy cutoff which is the flow parameter. We 
truncate the hierarchy by only considering the one- and 
two-particle vertices. The two-particle vertex is projected 
onto the Fermi points and parametrized by a static ef- 
fective interaction with local and nearest-neighbor parts. 
This implies a frequency independent S^. A detailed ac- 
count of our method was given in Refs.iniandEl Using 
the Dyson equation an approximate expression for Qi^n 
is obtained from taken at the end of the fRG flow at 
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FIG. 1: Upper panel: Conductance as a function of Vg for 
the SIAM at different U /A. Lower panel: Average number of 
electrons on the dot. 



A = 0. Generically the order N coupled differential equa- 
tions can only be integrated numerically. For a variety of 
transport problems through inhomogeneous LLs it was 
shown earlier that our method leads to reliable results 
for weak to intermediate interactions. Tl', T^, T?, 

Case (a) noninteracting leads, interacting dot: 
We focus on No = 1 for which our model reduces to the 
SIAM and here (for = 1) approximate the two-particle 
vertex by the bare interaction. Without LL leads, v only 
shifts the position of the resonance. We chose v = 1/2 
for which G(Vg) is symmetric around 0. The set of flow 
equations reduces to a single one for the effective onsite 
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with the initial condition V — Vg and the hybridiza- 
tion A — 27ri'2p, where p denotes the spectral weight at 
the end of the leads. As usual we here assume an energy 
independent p (infinite band width limit). Q Note that 
in this case U can be taken as the unit of energy. The 
upper panel of Fig.Qshows G{Vg) for different U/A. For 
[/ ^ A we recover the plateau-like resonance of unitary 
height. Also for asymmetric dot-lead couplings we re- 
produce the exact height of the plateau. The occupation 
of the dot can be computed from the Green function and 
is shown in the lower panel. In the plateau region it turns 
out to be close to 1 while it sharply raises/drops to 2/0 
to the left/right of the plateau. Our dot self-energy is 
frequency independent which leads to a Lorentzian dot 
spectral function of width 2 A and height l/(7rA) centered 
around V — V^^'^. This implies that the spectral weight 
at p, and thus GiVg) is determined by V.]^ The solution 
of the differential equation (j3Jl at A = is obtained in 
implicit form 



vJijv) - SJq{v) _ Jaivg) 
vYi{v) - SYo{v) ~ Yo{vg) 



(4) 



with V — Vtt/U, Vg — Vgir/U, S — Att/U, and Bessel 
functions J„, Yn- For \Vg\ < Vc, with Vc = VcTt/U being 
the first zero of Jq, i.e. Vc = 0.7655 U, this equation has a 
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FIG. 2: Upper panel: GiVg) for a noninteracting dot with 
LL leads at different A^. The parameters are No = 1, t' = 
0.1 1, n = 3/4, U — t, and U' = 0.65 1 (small two-particle 
backscattering) . Lower panel: Scahng of the resonance width. 
Filled circles: the same parameters as in the upper panel. 
Open circles: Nd = 1, i' = O.lt, n = 1/2, U ^ t, and U' = 
(large two-particle backscattering). 



solution with a small \V\. For U ^ A the crossover to a 
solution with \ V\ being of order U (for \ Vg\ > Vc) is fairly 
sharp. Expanding both sides of Eq. for small and 
\vg\ gives V — exp [— [//(ttA)]. The exponential pin- 
ning of the spectral weight at = for small \ Vg \ and the 
sharp crossover to a 1/ of order U when | Vg | > Vc leads to 
the observed resonance line shape. For U ^ A the width 
of the plateau is 2Vc — 1.531J7, which is larger than the 
width U found with NRG. 3j In the present context this 
difference is a minor issue. We expect that the agreement 
can be improved using a more elaborate fRG truncation 
scheme. |l6| It is remarkable that our technically fairly 
simple approximation reproduces the pinning of spectral 
weight, which is an important feature of Kondo physics. 
This has to be contrasted with other simple approxima- 
tions as e.g. low order perturbation theory or the self- 
consistent Hartree approximation that do not give the 
correct line shape of G{Vg). 

Case (b) LL leads, noninteracting dot: Using 
the fRG based approximation, we have earlier studied 
tunneling through a quantum dot embedded in a spinless 
LL. We showed that for ^ oo the resonance width w 
vanishes following a power law with a Kp dependent ex- 
ponent. Also the case of asymmetric dot- lead coupling 
was investigated. ^3 Including the spin degree of freedom 
the physics becomes more complex due to the possibil- 
ity of two-particle scattering of electrons with opposite 
spin at opposite Fermi points. This limits power laws 
to exponentially large length scales. To clearly observe 
LL behavior at experimentally accessible scales one has 
to consider a situation in which this backscattering pro- 
cess is small. In our model for fixed n and U this can 
be achieved by fine tuning t/'.^l^j In the upper panel of 
Fig. 121 we show the N dependence of G{Vg) for a sin- 
gle site dot computed for a small backscattering ampli- 
tude. At resonance voltage VJ the conductance is 2e^//i 



independent of N. In the lower panel the extracted w 
(filled circlea) is shown as a function of N. It follows the 
power law iV(^p-i)/2,0 with an fRG approximation to 
the LL parameter Kj, that is correct to leading order in 
the interaction. d [3 III El E.g. for n = U ^ t, 
and U' ~ 0.65 1 we find Kp = 0.76 in excellent agreement 
with the numerical result Kp ~ 0.7490.^ Off resonance 
G asymptotically vanishes oc N^~^f For Vg close 

to Vg such that 1 — G/{2e^/h) <^ 1, the deviation from 
the unitary limit scales as N^^^p, characteristic for a 
weak single impurity. Further increasing N this differ- 
ence increases and the behavior eventually crosses over 
to the off-resonance power-law suppression of G discussed 
above. Due to an exponentially large crossover scale, 
even for the very large N accessible with our method 
the complete crossover from one to the other power law 
cannot be shown for a single fixed Vg but follows from 
one-parameter scaling. |7t ilL] At sizeable backscattering 
the off-resonance conductance and thus the width w first 
slightly increase for increasing N - becoming larger than 
the noninteracting width - as shown by the open circles 
in the lower panel of Fig. |21 For larger N both quanti- 
ties start to decrease and eventually go to zero for expo- 
nentially large N. The backscattering process scales to 
zero (only) logarithmically in the low-energy limit and 
power-law behavior cannot be observed even for fairly 
long LL leads. (1^ An upper bound of the length of one- 
dimensional wires realized in experiments is of the order 
of /im, roughly corresponding to 10'' lattice sites, [l^ 

Case (c) — LL leads, interacting dot: We here fo- 
cus on the case in which the interactions on the dot and 
in the LL leads are taken to be equal. In the upper panel 
of Fig. 131 0(1/^) is shown for a parameter set with sizeable 
two-particle backscattering and LL lead length N = 10* 
typical for experiments. For interactions large compared 
to the hybridization we find the broad plateau-like res- 
onances induced by the Kondo effect, at least for finite 
LL leads. The same holds for other Nd, in particular 
for Nd — 1. The width of the plateaus is proportional 
to the local component of the effective interaction at the 
end of the fRG ffow and to 1/Nd- Here we are inter- 
ested in the interplay of Kondo and LL physics and thus 
focus on tunnel barriers with small transmission. In the 
plateau regions the number of electrons on the dot (lower 
panel) is odd while it is even for gate voltages where G 
is small. The upper panel of Fig. 0] shows the N de- 
pendence of G{Vg) computed for the same parameters 
as in the upper panel of Fig. El (small backscattering), 
but including the interaction on the dot. Note the dif- 
ferent X-axis scales of Figs. (21 and 01 In Fig. 01 the dif- 
ferences between the curves for different N are barely 
visible, in particular the changes of the resonance width 
are marginal. For parameters with sizeable backscatter- 
ing, as in Fig. O the difference between curves computed 
for different N are even smaller. To analyze the iV depen- 
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FIG. 3: Upper panel: G{Vg) 
leads. The parameters are: 
N = 10*, Nd =6,t' = O.lt. 
of electrons on the dot. 



for an interacting dot with LL 
n = 1/2, U = t, U' = 0.5t, 
Lower panel: Average number 
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FIG. 4: Upper panel: G{Vg) for an interacting dot with LL 
leads at different A'^. The parameters are as in the upper panel 
of Fig. 121 but with interaction on the dot. Lower panel: Scal- 
ing of G/{2e^/h) at Vg = outside the plateau (circles) and 
of 1 — G/{2e^/h) on the plateau at Vg — — 0.685 1 (squares). 



dence at small backscattering in more detail in the lower 
panel of Fig. 01 we show the scaling of G for a gate volt- 
age outside the plateau (circles) and of 2e^//i — G for a 
gate voltage on the plateau- like resonance (squares). For 
Vg outside the plateau G follows a power law with the 



exponent 1 — K~ and G vanishes for N ^ oo. Within 
every plateau we find a Vg at which G = 2e^//i indepen- 
dent of N. For Vg ^ Vg , still within the plateau, the 
deviation of G from the unitary limit scales as N^'^p, 
i.e. with the weak single impurity exponent. This shows 
that any deviation from Vg acts as an impurity. By anal- 
ogy with the single impurity behavior ^ll^ we conclude 
that in the asymptotic low-energy limit the impurity will 
effectively grow and for — > oo the plateaus will vanish. 
For infinitely long LL leads the resonances are infinitely 
sharp even in the presence of Kondo physics. However, 
for t' <^ U the plateaus at finite TV are well developed 
and the length scale on which they start to deteriorate is 
extremely large. 

Also for asymmetric dot-lead couplings we find (al- 
most) plateau-like resonances. To discuss this in more 
detail we focus on typical parameters with N « 10^ 
and an asymmetry A^/Aji « 2. Then the width is al- 
most unaffected by the asymmetry. For the interaction 



and filling as in Fig. 01 (small backscattering) the height 
within the plateaus varies by a few percent (with maxima 
at the left and right boundaries) and has average value 
w 0.85 (2e^//i). For sizeable backscattering the difference 
to the symmetric case is even smaller. With increasing 
A^ the variation of the conductance on the plateaus in- 
creases and the average value decreases. We expect that 
for N oo the resonances disappears. 

In summary, we studied how the linear conductance 
through a quantum dot is modified in the presence of 
both Kondo physics and LL leads. Using an approxi- 
mate method that is based on the fRG we investigated 
the dependence of G on the gate voltage as well as the 
length A^ of the LL leads. We found that for all experi- 
mentally accessible length scales and for typical left-right 
asymmetries of the dot-lead hybridizations the plateau- 
like resonances characteristic for Kondo physics will also 
be present if the leads are LLs. The plateaus are more 
pronounced if the two-particle backscattering is sizeable, 
although they disappear for A^ — > oo. 
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